Simulação numérica da gaseificação de biomassa em leito fixo

O objetivo deste trabalho foi avaliar a potencialidade das sementes de açaí rejeitadas após beneficiamento em Belém – Pará, a fim de ser usada como fonte de biomassa para a produção de energia por gaseificação. Resultados da análise elementar obtidas no LABCAT – UNB em base seca foram obtidas com percentuais de Carbono, 46,06 %, Hidrogênio 6,01 % Oxigênio 43,37 %, Nitrogênio 4,33 % e Enxofre 0,08%. Tais valores entre outros parâmetros obtidos em estudo experimental realizados por Santos 2011 e Cruz, 2010, foram utilizados para realizar uma simulação numérica do processo de gaseificação por meio de um código FORTRAN zero dimensional desenvolvido especificamente para esta finalidade. Variações da influência da razão de equivalência e teor de umidade nas concentrações de CO, H2 e CH4 foram comparados a resultados experimentais. Um modelo transiente unidimensional também foi analisado. Utilizando a aplicação de um mecanismo de pirólise de dois passos onde o alcatrão primário (ou oxigenados) são inicialmente formados e então craqueados em alcatrão secundário (ou hidrocarbonetos) e outros gases combustíveis. Assim se revela a complexidade dos processos físicos e químicos que ocorrem no reator pelo uso do balanço de energia e massa e com informações das taxas de reação das reções químicas e os processos de transporte físicos. O modelo computacional é capaz de prever o perfil da composição de gases, temperatura, alcatrão primário e secundário ao longo do eixo longitudinal do gaseificador. Foram obtidos resultados para simulações realizadas nos dois modelos e posteriormente analisado tal comportamento em comparação ao resultado experimental obtidos por Santos (2011).

A manutenção do equilíbrio econômico de contratos de obras públicas empreitadas: parâmetros para o estudo de riscos vinculados a atrasos no encaixe do preço

Dentre os diversos segmentos do Setor da Construção Civil, um dos principais é o de obras por empreitada global, que são aquelas contratadas por preço certo e total. As obras públicas, que têm como finalidade atender a uma utilidade pública e cujos contratos são regidos pela Lei 8.666/93, são geralmente contratadas sob esse regime de execução. O artigo 58 desta Lei determina que deva ser mantido o equilíbrio econômicofinanceiro dos contratos celebrados com o Poder Público. A fundamentação teórica trata dos conceitos de obras públicas, dos principais aspectos da Lei 8.666/93, de análise econômico-financeira e de análise de riscos. Este trabalho apresenta também uma sistemática e uma simulação numérica para formação de preços para obras empreitadas de Construção Civil. O objetivo geral deste trabalho é a criação de parâmetros para o estabelecimento de margens de proteção para a dilatação do prazo e respectivos encaixes de preços, visando, como determina a Lei 8.666/93, a manutenção do equilíbrio original. Com base na análise dos resultados de 63 obras, verifica-se a manutenção do equilíbrio econômico-financeiro, comparando-se a taxa de retorno restrita esperada pelas empresas no momento das licitações com a realmente alcançada na operação. A partir dos desvios detectados nas taxas de retorno nesses dois momentos, cria-se, pelo processo de simulação e análise de sensibilidade, parâmetros para cobertura dos referidos riscos. No momento da licitação, a taxa de retorno restrita média das obras analisadas é de 9,45% ao mês. No segundo momento, essa taxa cai para em média para 5,16% ao mês. Isto representa, em média, uma variação de 45,39%. A taxa média de atratividade das obras públicas analisadas está dentro do intervalo de 6% a 13% ao mês. Como parâmetro de cobertura de riscos quanto aos desvios de prazo e encaixe do preço, tem-se o intervalo de 1,5% a 11,5% do preço. Portanto, entre as obras analisadas, não existe o equilíbrio econômico-financeiro.

Estudo dos potenciais impactos das mudanças climáticas e de alterações na cobertura vegetal nos recursos hídricos na Região Central da Amazônia

Este trabalho busca avaliar e quantificar os impactos na disponibilidade hídrica local na região da Amazônia Central, decorrentes de possíveis mudanças climáticas e modificações na cobertura vegetal, por meio de um experimento de simulação numérica com o modelo de biosfera Common Land Model (CLM), no modo “off line”. Os resultados de 9 modelos climáticos acoplados Oceano-Atmosfera para três cenários de alterações climáticas do Painel Intergovernamental sobre Mudanças Climáticas (IPCC-AR4) foram utilizados para compor a base das forçantes climáticas do modelo de biosfera CLM. Em relação ao uso da terra, foram utilizados cenários de dinâmica de desmatamento para a região no caso “business as usual” previstos para cada ano do período de 2001 a 2050. A área de estudo compreende o domínio da drenagem da Bacia do Rio Cuieiras na Amazônia Central. A partir dos resultados dos modelos analisaram-se, para cada conjunto de simulações, as incertezas das projeções em relação à precipitação e a temperatura e o impacto no ciclo hidrológico terrestre, considerando a variabilidade entre os modelos e os cenários de emissões de CO2; bem como as alterações nas componentes do balanço de água e energia a superfície associada a variações progressivas na cobertura de floresta e sua substituição por pastagem. Os resultados indicam que diante de um cenário de mudanças climáticas que resultem em uma diminuição (aumento) persistente na precipitação média anual, tanto o escoamento quanto o armazenamento de água no solo serão diretamente afetados. Em relação às alterações na cobertura vegetal tanto as componentes dos balanços de água e energia foram significativamente afetados pela substituição da floresta por pastagem, apresentando redução na evapotranspiração e aumentos no armazenamento de água no solo e do escoamento total.

Análise computacional dos campos de velocidade e temperatura do ambiente interno da usina termelétrica Santana - Amapá

A simulação numérica do escoamento de ar em ambientes internos é na atualidade o método mais apropriado para análise de conforto térmico em ambientes internos. O escoamento de ar nesses ambientes configura-se como um escoamento complexo, pois, em regra geral, é uma combinação de escoamentos cisalhantes livres (jatos) e cisalhantes de parede, além disso, esses escoamentos são governados por forças de inércia e forças de empuxo, caracterizando-o como de convecção mista. A combinação desses mecanismos cria um escoamento com características complexas, como zonas de recirculação, vórtices, descolamento e recolamento de camada-limite dentre outras. Portanto, a precisão da solução estará diretamente ligada, principalmente, na habilidade do modelo de turbulência adotado de reproduzir as características turbulentas do escoamento de ar e da transferência térmica. O objetivo principal do presente trabalho foi a simulação computacional do ambiente térmico interno do galpão que abriga os geradores e motores Wärtzilä da Usina Termelétrica Santana no estado do Amapá. A formulação matemática baseada na solução das equações gerais de conservação inclui uma análise dos principais modelos de turbulência aplicados ao escoamento de ar em ambientes internos, assim como os processos de transferência de calor associados. Na modelagem numérica o método de volumes finitos é usado na discretização das equações de conservação, através do código comercial Fluent-Airpak, que foi usado nas simulações computacionais para a análise dos campos de velocidade e temperatura do ar. A utilização correta do programa computacional foi testada e validada para o problema através da simulação precisa de casos retirados da literatura. Os resultados numéricos foram comparados a dados obtidos de medições experimentais realizados no galpão e apresentou boa concordância, considerando a complexidade do problema simulado, o objetivo da simulação em face da diminuição da temperatura no interior do galpão e, também, em função das limitações encontradas quando da tomada das medições experimentais. Além disso, foram feitas simulações de estratégias de melhoria do ambiente térmico da Usina, baseadas na realidade levantada e nos resultados da simulação numérica. Finalmente, foram realizadas simulações do protótipo de solução proposto para a diminuição da temperatura interna do galpão o que possibilitará um aumento, na faixa de 20 a 30%, do tempo de permanência no interior do galpão.

Projeto e construção de uma câmara reverberante em escala reduzida para a caracterização acústica de materiais absorventes

Na área da Acústica, no que diz respeito às últimas décadas, tem crescido substancialmente a demanda por informações sobre a capacidade de absorção dos materiais acústicos, bem como a utilização de métodos de modelagem que permitam inferir o comportamento acústico dos sistemas reais, os quais demandam uma infraestrutura de porte, o que representa um alto investimento. Assim, visando atender a estas demandas, possibilitar o desenvolvimento de pesquisas sobre as características de absorção sonora de materiais tipicamente regionais e o desenvolvimento de aulas práticas no âmbito da acústica, tanto no curso de graduação quanto na pós-graduação, este trabalho apresenta todas as informações pertinentes ao projeto e construção de uma mini-câmara reverberante, na escala de 1: 6 de um volume de 200 m³. Além disso, é apresentada uma metodologia de qualificação e adequação do campo acústico interno, tendo por base as orientações contidas na norma ISO-354. Por outro lado, são apresentados e discutidos os resultados obtidos, a partir de ensaios realizados nesta mini-câmara, caracterizando e validando os painéis de fibra de coco, desenvolvidos pelo Grupo de Vibração e Acústica, da Universidade Federal do Pará, bem como uma simulação numérica por Raios Acústicos, tendo por base a cavidade acústica da minicâmara, para permitir a validação do modelo numérico a partir dos dados experimentais.

Modelagem e filtragem de distorções estáticas em levantamentos de dipolos contínuos

As distorções estáticas são um grave problema que afeta o método magnetotelúrico e muitas tentativas têm sido feitas para eliminar ou minimizar os seus efeitos. Este trabalho trata de uma técnica nova que trata deste problema, o método EMAP, que é uma adaptação do método magnetotelúrico na qual as medidas de campo elétrico são feitas em uma linha contínua de dipolos conectados entre si, e os dados coletados desta maneira são tratados com um filtro espacial passa-baixa dependente da frequência. Este trabalho é composto de duas partes principais, a simulação numérica dos dados contaminados com as distorções estáticas e a filtragem destes dados com o filtro espacial passa-baixa. Na primeira parte, aplicamos o método dos elementos finitos para simular a resposta dos dipolos elétricos, os quais fazem as medidas do campo elétrico. Na segunda parte aplicamos a janela de Hanning como filtro passa-baixa apropriado para tratar os dados.

Simulação de perfis nucleares de poço em formações complexas

A identificação e descrição dos caracteres litológicos de uma formação são indispensáveis à avaliação de formações complexas. Com este objetivo, tem sido sistematicamente usada a combinação de ferramentas nucleares em poços não-revestidos. Os perfis resultantes podem ser considerados como a interação entre duas fases distintas: • Fase de transporte da radiação desde a fonte até um ou mais detectores, através da formação. • Fase de detecção, que consiste na coleção da radiação, sua transformação em pulsos de corrente e, finalmente, na distribuição espectral destes pulsos. Visto que a presença do detector não afeta fortemente o resultado do transporte da radiação, cada fase pode ser simulada independentemente uma da outra, o que permite introduzir um novo tipo de modelamento que desacopla as duas fases. Neste trabalho, a resposta final é simulada combinando soluções numéricas do transporte com uma biblioteca de funções resposta do detector, para diferentes energias incidentes e para cada arranjo específico de fontes e detectores. O transporte da radiação é calculado através do algoritmo de elementos finitos (FEM), na forma de fluxo escalar 2½-D, proveniente da solução numérica da aproximação de difusão para multigrupos da equação de transporte de Boltzmann, no espaço de fase, dita aproximação P₁, onde a variável direção é expandida em termos dos polinômios ortogonais de Legendre. Isto determina a redução da dimensionalidade do problema, tornando-o mais compatível com o algoritmo FEM, onde o fluxo dependa exclusivamente da variável espacial e das propriedades físicas da formação. A função resposta do detector NaI(Tl) é obtida independentemente pelo método Monte Carlo (MC) em que a reconstrução da vida de uma partícula dentro do cristal cintilador é feita simulando, interação por interação, a posição, direção e energia das diferentes partículas, com a ajuda de números aleatórios aos quais estão associados leis de probabilidades adequadas. Os possíveis tipos de interação (Rayleigh, Efeito fotoelétrico, Compton e Produção de pares) são determinados similarmente. Completa-se a simulação quando as funções resposta do detector são convolvidas com o fluxo escalar, produzindo como resposta final, o espectro de altura de pulso do sistema modelado. Neste espectro serão selecionados conjuntos de canais denominados janelas de detecção. As taxas de contagens em cada janela apresentam dependências diferenciadas sobre a densidade eletrônica e a fitologia. Isto permite utilizar a combinação dessas janelas na determinação da densidade e do fator de absorção fotoelétrico das formações. De acordo com a metodologia desenvolvida, os perfis, tanto em modelos de camadas espessas quanto finas, puderam ser simulados. O desempenho do método foi testado em formações complexas, principalmente naquelas em que a presença de minerais de argila, feldspato e mica, produziram efeitos consideráveis capazes de perturbar a resposta final das ferramentas. Os resultados mostraram que as formações com densidade entre 1.8 e 4.0 g/cm³ e fatores de absorção fotoelétrico no intervalo de 1.5 a 5 barns/e^-, tiveram seus caracteres físicos e litológicos perfeitamente identificados. As concentrações de Potássio, Urânio e Tório, puderam ser obtidas com a introdução de um novo sistema de calibração, capaz de corrigir os efeitos devidos à influência de altas variâncias e de correlações negativas, observadas principalmente no cálculo das concentrações em massa de Urânio e Potássio. Na simulação da resposta da sonda CNL, utilizando o algoritmo de regressão polinomial de Tittle, foi verificado que, devido à resolução vertical limitada por ela apresentada, as camadas com espessuras inferiores ao espaçamento fonte - detector mais distante tiveram os valores de porosidade aparente medidos erroneamente. Isto deve-se ao fato do algoritmo de Tittle aplicar-se exclusivamente a camadas espessas. Em virtude desse erro, foi desenvolvido um método que leva em conta um fator de contribuição determinado pela área relativa de cada camada dentro da zona de máxima informação. Assim, a porosidade de cada ponto em subsuperfície pôde ser determinada convolvendo estes fatores com os índices de porosidade locais, porém supondo cada camada suficientemente espessa a fim de adequar-se ao algoritmo de Tittle. Por fim, as limitações adicionais impostas pela presença de minerais perturbadores, foram resolvidas supondo a formação como que composta por um mineral base totalmente saturada com água, sendo os componentes restantes considerados perturbações sobre este caso base. Estes resultados permitem calcular perfis sintéticos de poço, que poderão ser utilizados em esquemas de inversão com o objetivo de obter uma avaliação quantitativa mais detalhada de formações complexas.

Modelagem e simulação do processo de separação por membrana difusa polimérica para a extração de metais

O processo de extração de metais em membrana polimérica é analisado empregando o método das Linhas-Gear nas equações diferenciais resultantes da modelagem matemática do problema físico. O estudo tem a finalidade de averiguar o comportamento das curvas de extração de metal extraído e a influência dos parâmetros físicos da membrana no processo de extração por membrana, e também analisar, a partir de um estudo paramétrico, a influência dos parâmetros mais importantes do processo, verificando-se assim se o processo de extração se encontra sob controle cinético ou difusivo. A modelagem do problema é realizada assumindo-se a extração dos metais cádmio, Cd (II), ouro, Au (lII), e Paládio, Pd (11) sobre a membrana difusiva polimérica Aliquat 336/PVC. Os resultados são então comparados com dados obtidos em trabalhos previamente reportados na literatura, concluindo-se que os mesmos são satisfatórios.

Desenvolvimento de Interface gráfica como suporte para soluções numéricas das equações de Maxwell em coordenadas gerais – 3D

Neste trabalho, é implementada uma interface gráfica de usuários (GUI) usando a ferramenta Qt da Nokia (versão 3.0). A interface visa simplificar a criação de cenários para a realização de simulações paralelas usando a técnica numérica Local Nonorthogonal Finite Difference Time-Domain (LN-FDTD), aplicada para solucionar as equações de Maxwell. O simulador foi desenvolvido usando a linguagem de programação C e paralelizado utilizando threads. Para isto, a biblioteca pthread foi empregada. A visualização 3D do cenário a ser simulado (e da malha) é realizada por um programa especialmente desenvolvido que utiliza a biblioteca OpenGL. Para melhorar o desenvolvimento e alcançar os objetivos do projeto computacional, foram utilizados conceitos da Engenharia de Software, tais como o modelo de processo de software por prototipagem. Ao privar o usuário de interagir diretamente com o código-fonte da simulação, a probabilidade de ocorrência de erros humanos durante o processo de construção de cenários é minimizada. Para demonstrar o funcionamento da ferramenta desenvolvida, foi realizado um estudo relativo ao efeito de flechas em linhas de baixa tensão nas tensões transitórias induzidas nas mesmas por descargas atmosféricas. As tensões induzidas nas tomadas da edificação também são estudadas.

Simulação de sistemas de simples estágios de refrigeração por compressão de vapor

Na atualidade, o estudo do desempenho térmico de um sistema de refrigeração por compressão de vapor representa uma ferramenta importante no auxílio do desenvolvimento de novos produtos ou melhoria dos já existentes. Um modelo de simulação em regime permanente foi elaborado para avaliar o desempenho do sistema frigorífico. O sistema estudado inclui uma Central de Ar Condicionado, modelo PA HILTON, constituída de um compressor alternativo do tipo semi-hermético, evaporador e condensador compacto de tubos e aletas e uma válvula de expansão termostática. O modelo do condensador considera três regiões distintas de troca de calor as quais são respectivamente a região de dessuperaquecimento, condensação e subresfriamento. Para a modelagem do evaporador, foram consideradas as regiões de evaporação e superaquecimento. No modelo de simulação foram utilizadas correlações adequadas para a estimativa dos coeficientes de transferência de calor e perda de pressão para cada região do evaporador e condensador. Não foram consideradas a transferência de calor e queda de pressão nas linhas de conexão entre os componentes. A solução do sistema de equações não lineares resultantes da modelagem matemática dos componentes do sistema simulado foi obtida utilizando-se o método das substituições sucessivas com o emprego do software Engineenng Equation Solver . Os resultados obtidos pelo modelo de simulação apresentaram erros inferiores a 9% em relação aos valores experimentais.

Simulação de transporte eletrônico em dispositivos unimoleculares baseados em indicadores de pH

Nas últimas décadas, diversos pesquisadores têm tentado empregar moléculas em dispositivos eletrônicos de nanoescala. Por este motivo, diferentes parâmetros eletro/ópticos, que regem o transporte eletrônico em moléculas orgânicas, precisam ser analisados. Neste trabalho foi desenvolvido um estudo de transporte de carga para o composto Vermelho de Propila, popularmente utilizado como indicador de pH. A motivação para estudá-lo resulta de sua estrutura constituída por subunidades doadora-aceitadora, acopladas via grupo azo (N=N), uma característica bem conhecida em retificadores moleculares. A metodologia utilizada para tratar o sistema em equilíbrio é baseada em métodos de Mecânica Molecular e Hartree-Fock. Sendo que, para simular o sistema em não-equilíbrio, foi empregado o formalismo de Landauer-Büttiker. Através desses métodos, as curvas características do sistema molecular foram traçadas e comparadas. O resultado da comparação permitiu explicar os fenômenos que regem o transporte eletrônico na nanoestrutura. Além disso, foram analisados os efeitos de contatos metálicos, ligados a molécula na presença de campo elétrico externo.

Metodologia de discretização espacial para o método radial de interpolação por pontos (RPIM) aplicada para solução numérica das equações de Maxwell

Neste trabalho, foi desenvolvido e implementado um método de discretização espacial baseado na lei de Coulomb para geração de pontos que possam ser usados em métodos meshless para solução das equações de Maxwell. Tal método aplica a lei de Coulomb para gerar o equilíbrio espacial necessário para gerar alta qualidade de discretização espacial para um domínio de análise. Este método é denominado aqui de CLDM (Coulomb Law Discretization Method ) e é aplicado a problemas bidimensionais. Utiliza-se o método RPIM (Radial Point Interpolation Method) com truncagem por UPML (Uniaxial Perfectlly Matched Layers) para solução das equações de Maxwell no domínio do tempo (modo TMz).

Desenvolvimento de framework para análise e simulação dinâmica de sistemas elétricos de potência

Esta tese apresenta uma metodologia flexível orientada a objetos (OO) para a aplicação no projeto e implementação de sistemas de software utilizados na realização de estudos dinâmicos de sistemas elétricos de grande porte. A metodologia OO proposta objetiva tornar mais simples o desenvolvimento, a atualização e a manutenção de complexos sistemas de software para estudos de transitórios eletromecânicos em sistemas elétricos de potência. Os requisitos de usuário são mapeados para um conjunto de classes básicas, as quais são usadas para efetuar a modelagem de dispositivos dinâmicos tais como geradores elétricos. Para avaliação da metodologia foram realizados dois estudos de casos. No primeiro estudo caso o Framework foi aplicado na simulação das unidades geradoras da Usina Hidrelétrica de Tucuruí. Os resultados da simulação foram comparados com medições obtidos em ensaios no campo e mostrou a boa performance do Framework na reprodução dos fenômenos eletromecânicos desta usina de grande porte. No segundo estudo de caso, por outro lado, o Framework foi aplicado na modelagem de um sistema de geração fotovoltaico (PV) com seu sistema de Rastreamento da Potência Máxima (MPPT). O controle MPPT foi implementado usando técnicas digitais. Os resultados das simulações demonstram a performance do Framework na modelagem do sistema de controle de corrente, assim como no controle MPPT, dos sistemas de geração PV.

Simulação computacional no auxílio à tomada de decisão nos processos industriais: utilizando a ferramenta Tecnomatix Plant Simulation 9.0

Atualmente existem inúmeros processos produtivos automatizados, os quais vêm se tornando cada vez mais complexos em função das necessidades do mundo moderno e, portanto, demandam nas fases de projetos e de implementação ferramentas de engenharia cada vez mais poderosas para modelá-los e analisá-los da maneira mais eficiente possível. Nesse ambiente de crescente pressão por resultados positivos, racionalização e aprimoramento de recursos internos, que a ferramenta computacional Tecnomatix Plant Simulation 9.0, surge como um caminho para a obtenção de competitividade produtiva. Ressaltando ainda que, o estudo proposto é de grande relevância para os profissionais da gestão produtiva, os quais almejam resultados que minimizem custos e maximizem lucros.

Identificação de correntes de inrush na proteção diferencial de transformadores de potência através do gradiente da corrente diferencial e de mapas auto-organizáveis

A principal dificuldade encontrada na proteção diferencial de transformadores de potência é a correta distinção entre as correntes de inrush e as correntes de faltas internas. Tradicionalmente os relés diferenciais executam esta tarefa utilizando a técnica de restrição por harmônicos baseada na premissa de que as correntes de inrush possuem alta concentração de componentes harmônicas de segunda ordem, contudo essa técnica nem sempre é eficaz. O presente trabalho tem como objetivo apresentar a proposta de duas novas metodologias capazes de realizar a identificação e distinção entre as correntes de inrush das correntes de faltas internas na proteção diferencial de transformadores de potência através de metodologias que não dependem do conteúdo de harmônicos do sinal da corrente diferencial. A primeira metodologia proposta, denominada de método do gradiente da corrente diferencial, é baseada no comportamento do vetor gradiente, obtido através da diferenciação numérica do sinal da corrente diferencial. O critério de distinção utilizado é baseado no desvio padrão do ângulo do vetor gradiente que apresenta comportamento diferenciado para correntes de inrush e correntes de curto-circuito. A segunda metodologia proposta é baseada na capacidade de reconhecimento e classificação de padrões das redes neurais de Mapeamento Auto-organizável de Kohonen. Como padrão de entrada e de treinamento da rede neural é utilizado um vetor contendo quatro níveis do espectro do desvio padrão do ângulo do vetor gradiente da corrente diferencial nas três fases do transformador de potência. A eficácia dos métodos foi testada através da simulação de diversas situações de faltas internas e correntes de inrush, incluindo situações de “Sympathetic Inrush”, em um transformador de potência usando o software EMTP/ATP e através da implementação do algoritmo em MATLAB®, apresentando resultados altamente promissores.